Generic f(R) theories and classicality of their scalarons 
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In this Letter, we study quantum stability bound on the mass of scalaron in generic theories of 
f(R) gravity. We show that in these scenarios, the scalaron mass increases faster with local density 
of the environment than one-loop quantum correction to it thereby leading to violation of quantum 
bound on the chameleon mass. The introduction of quadratic curvature corrections in the action 
are shown to stabilize the model. 



I. INTRODUCTION 

The late time cosmic acceleration 0, Q has recently 
been accepted as one of the fundamental phenomena of 
nature whose underlying cause remains yet to be un- 
folded. The standard lore preaches that the late time 
acceleration is caused by the presence of a cosmic fluid 
with large negative pressure; the cosmological constant 
A 0] presents a distinguished example of such a fluid. 

As an alternative to cosmological constant, a variety 
of scalar field models were investigated with a hope 
to alleviate the fine tuning and coincidence problems 
associated with the model without assigning a funda- 
mental reason to switch off A. Unfortunately, the scalar 
field dark energy models are not entirely problem free, 
assumptions about model parameters/tuning are tacitly 
made in these models. 

There is an alternative school of thought in cosmology 
which advocates the need for paradigm shift and believes 
that cosmic acceleration results from large scale modifi- 
cation of gravity Q. Such a proposal sounds healthy as 
general theory of relativity, which passes the solar test 
with great precision, is often extrapolated to large scales 
where it was never verified directly. 

One of the schemes of large scale modification based 
upon phcnomcnological considerations is provided by 
f(R) theories of gravity [!, [||. These theories essen- 
tially contains an additional scalar degree of freedom 
apart from graviton. Indeed, f(R) theories are con- 
formally equivalent to Einstein theory plus a canonical 
scalar degree of freedom dubbed scalaron whose poten- 
tial is uniquely constructed from Ricci scalar. 

It is interesting to note that f(R) gravity is Os- 
trogradski ghost free despite the equations of motion 
being of fourth order as there are enough number of 
constraints to protect the theory. 



*On leave from Centre for theoretical physics, Jamia Millia islamia, 
New Delhi 



However, it should further be ensured that the gravi- 
ton and the scalaron are well behaved which imposes re- 
strictions on the functional form of f(R). Namely, the 
generic f(R) theories should satisfy, f'(R),f"(R) > 
in order to avoid the ghost and tachyonic modes. Sec- 
ondly, these theories should reduce to ACDM in the high 
density regime in order to comply with the local gravity 
constraints. The class of models proposed by Hu-Sawicki 
and Starobinsky (HSS) 0,1a] (see also Ref.Q on the same 
theme) satisfy the said requirements and are of great in- 
terest in f(R) theories. 

In this scenario, the scalar degree of freedom is non- 
minimally coupled to matter in Einstein frame and hence 
it might conflict with the local physics which does not see 
a fifth force. Thus if all is to be well, the scalaron should 
acquire a heavy mass in local environment in order to 
suppress the fifth force and become light with mass of 
the order Hubble constant today to be relevant to late 
time cosmic acceleration which means that the scenario 
asks for a chameleon. 

In f(R) theories, the scalaron mass naturally acquires 
density dependence thereby allowing us to naturally 
implement the chameleon mechanism by appropriately 
choosing the form of f(R) giving rise to higher values 
of scalaron mass for larger values of density of the 
environment. The chameleon scenario, despite of its 
attractiveness, is plugged with several difficult problems: 
The scalaron mass might exceed the Planck mass by sev- 
eral orders of magnitudes in high density configurations 
such as neutron stars, the curvature singularity is easily 
accessible in the scenario and requires ugly fine tuning 
for its cu rellOl . Being inspired by Starobinsky 's original 
proposal pjj, the HSS model was extended by adding 
quadratic curvature correction [l2[ to address the said 
problems. 

The quadratic correction provides in a sense quantum 
correction to gravity sector which turns out to be im- 
portant in the scenario under consideration. It becomes 
equally important to investigate whether the quantum 
1-loop correction to scalaron potential remains small as 
density of the environment increases. 

In this Letter we shall study the quantum stability 
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bound for scalaron in Starobinsky f(R) gravity model. 
We also address the same issue in the framework of an 
extended scenario by incorporating the quadratic curva- 
ture corrections in the Starobinsky model. 



Chameleon field 

Let us consider the following action in the Einstein 
frame 



r M 2 1 
S = I d'x^—gy-^R - -(V^) 2 - V(4>) 



+5„ 



(1) 



The equation for the held </> which follows from the 
action can be written as 



(2) 



where {p, P) are energy density and the pressure in 
the Jordan frame. We consider this frame as our physi- 
cal frame in which the stress-energy tensor is conserved 
hence we assume that our quantities are independent of 
the scalar held <j>. We note that in the original paper [l4| , 
the authors defined a conserved density in the Einstein 
frame for a FLRW space-time. The definition that we 
adopt here gives a definition of the effective potential for 
any background (also in presence of pressure) and within 
this definition the effective mass of the chameleon field is 
the mass of the scalaron in f(R). It is however clear that 
because in most of the cases A ~ 1 the quantities in the 
two frames are very close. 

The eq. @ can be cast in the form 



n<p = 



dV f 



eff 



(3) 



L eff 



A' A" 

'4+4- 



15) 



To avoid a ghost instability, we require that V" /V > 

A A' ' 

In the simple case when A is given by, A = e l3c ^^ Mp , we 
have 



l eff 



V" 



-V 



(6) 



In what follows we shall consider the chameleon mecha- 
nism in f(R) theories of gravity where it occurs naturally. 

A. Chameleon theory and f(R) gravity 

Let considering f(R) action in the Jordan frame, 
Ml 



M f 

S =^f L J d i x^f(R)+S m [g^ i m i \. 
We next use a conformal transformation 



9nv -> f,R 9V = e 



-2l3(j>/M pl 



(7) 



(8) 



with j3 = — 1/V6, to transform the action to the Ein- 
stein frame 



where 



d X\J ~g 



(9) 



(10) 



where V eff = V+ ^§^A A . 

The existence of the chameleon mechanism depends on 
the form of the effective potential which in turn depends 
on the local density and pressure. When pressure is neg- 
ligible and density is large, the scalar field may acquire 
a large mass for a suitably chosen potential leading to 
suppression of the fifth force locally. The scalaron mass 
is defined as usual 



d 2 V e 



L eff 



eff 



(4) 



The scalar field is assumed to be settled in the mini- 
mum of the effective potential. It is therefore simple to 
recast the effective mass in the following form 



The hrst and second derivatives of the potential V((f>) 
are given by 



V.rhr. 



M pl 2/ - Rf. R 

Ve f 2 R ' 



(ii) 



3/, 



RR 



^ Rf.RR _ 4//.HJ? N | 



f,R 



One can see that effective potential belongs to 
Chameleon theory as {<p directly couples to matter), 



Veff = V{4>) + 



(13) 



where p and P are respectively the density and the pres- 
sure in Jordan frame and A = 1 / -J f t R. 
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It is easy to find that the minimum of the effective 
potential from cq. (|llll3[) . 



In |l5|, the authors considered the one- loop Coleman- 
Weinberg correction to the chameleon field potential. 



2/ - Rf, R = 



p-3P 
M 2 



(14) 



It is interesting to notice that the minimum of the po- 
tential is invariant under the addition to the action of a 
i? 2 -term. We shall use this aspect in the discussion to 
follow. 

Also one can rewrite the effective mass © with the 
help of (|11I12|) as 



l eff 



3f,RR 



1 - 



RR 



f,R 



(15) 



which corresponds to the mass of the scalaron in 
the Einstein frame. In fact in the Jordan frame, we 

f.R - Rf.RR} and therefore M = 



have M 2 = rr-r — 

OJ,RR 

\J f,R m eff = A~ 1 m e ff which is the standard factor 
which relates the mass of the field in Jordan to its counter 
part in Einstein frame. We should emphasis that we re- 
cover the mass of the scalaron because we consider the 
effective potential instead of the potential and also be- 
cause we consider the density p and the pressure in the 
Jordan frame as defined in (1131) . 

We often encounter local densities much larger than 
cosmological density p cr such as pi a b ~ 10 g/cm 3 and use 
the classical description for scalar degree of freedom in 
f(R) which assumes the quantum correction to scalaron 
potential to be small. Following Ref.[l5|, we shall now 
address the issue of quantum stability in generic theories 
of f(R) gravity. 



L eff 



64tt 2 



In 



°eff 



(</>) 



Mo 



(16) 



where fio is a cut off mass scale. It can be chosen 
equal to the mass of the field for a particular environ- 
ment (density) which would kill the quantum correction 
but the correction would reappear at other values of the 
density. 

At large values of density of interest to us or large mass 
of the field, we can set log term to unity 



" l eff 

64tt 2 ' 



(17) 



Since we expect quantum corrections to be small, 
we should have [l5[ a small modification of the shape 
of the potential V which implies AVi—i 0O p,^/V^ and 
AV\-ioop,<t><t>jy,^ should be small. 

Secondly, at the minimum of the effective potential 
(fT5]l, we have 2 



V' (0) + JL p e 4 W M P > = 

M„r 



(18) 



from which we obtain 



dp 



M p im 2 eff {4>) 



4W/M P) 



(19) 



which gives 



Quantum stability bound 

The effective potential defined in ([T3"|) depends on the 
energy density and therefore on the position of space- 
time. At the equilibrium, the field minimizes the poten- 
tial Veff.ff, = 0, hence the chameleon appears as a mas- 



"eff 



(j) /2).In Einstein frame, we have 



sive field (V e ff 

General Relativity with a scalar field. In scenario under 
consideration, the curvature scalar is small as seen later 
and the effects of the expansion are negligible. Hence the 
model is close to quantum field theory in flat space-time 
where we can neglect the effects of gravity 1 and quantize 
the scalar field sector in standard way. Also we should 
emphasize that we can always work in the Einstein frame, 
as long as the conformal transformation is not singular. 
Indeed this the case for the model studied as R ~ p/M 2 ^ 
to be demonstrated shortly. 



1 The gravity sector is the massless spin two particle without the 
scalaron. 



AV\-loop,cl> 



M 



/-2 



1 dm 6 



eff 



96tt 2 /3 2 p dp 



< 1 



Ml d 2 - 6 



eff 



96ir 2 /3 2 dp 2 

and after integration 



< 1 



48tt 2 /? 2 P 2 



m eff < 



1/6 



0.0073 



(20) 
(21) 



lOg cm 



-3 



1/3 



eV. 



(22) 



where, the constant of integration is fixed to zero; we 
can demand that the correction is zero for very low den- 
sities where the effective mass is zero. Let us briefly 
comment on the viability of Coleman- Weinberg one loop 



For simplicity we neglect the pressure 
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correction used here for chameleon. It corresponds to the 
quantum mechanically corrected potential 



V t , 



V eff ^) + -\ndet[d 2 



)] 



(23) 



where the first term represents classical part of the 
potential. The quantum correction is formally divergent 
and requires ultraviolet cut off. In case we use non- 
covariant scheme of regularization, say, Pauli-Willars 
regularization, with cut off M uv , the quantum correction 
is represented by three terms: (1)M„„, (2) m\^M 2 v and 
(3) ln(m^y). The first term can be absorbed in the 
definition of renormalization of cosmological constant, 
the third term represents the one-loop quantum correc- 
tion to be used in the discussion to follow. However, the 
second term is much larger than the second and would 
invalidate usage of the quantum bound based upon the 
third term only. 



bound obtained using the Coleman- Weinberg formula for 
effective potential. 

In f(R), (3 = — l/v6j which implies that 



m eff < 5.4 1(T 3 H 

\ lOg cm" 



1/3 



eV. 



(26) 



Eq. ([26|) provides an upper bound on the mass of the 
field. As we mentioned before, we shall be interested in 
the scrutiny of generic f(R) theories, namely HSS and 
would specialize to Starobinsky parametrization for con- 
venience. 



Parameters of Starobinsky model 

We are interested to study the quantum stability of 
Starobinsky f(R) gravity [1] 



It is well known that the term quadratic in cut off 
is specific to any regularization scheme which breaks 
Lorentz symmetry. In case of gauge theories, the reg- 
ularization scheme which does not respect the underly- 
ing symmetry of the theory leads to wrong results fl7l ]. 
Indeed, in the present context, the regularized value of 
the quantum correction using dimensional regularization 
gives rise to expression (|16[) without the dangerous term 
quadratic in cut off. It is interesting to note that we see 
similar features when we regularize the vacuum energy. 
In fact the correction can be understood as the zero-point 
energy density of the scalaron. It appears as a massive 
Klein- Gordon field and gives for the energy density of the 
vacuum 



f(R) = R- fiR c 1 - (1 + R 2 /R 



(27) 



During the de-Sitter phase, the solution is described by 
([M]) in an empty Universe. The curvature scalar (Rds) 
is solution of 



1 

M = o 



x(l + x 2 ) 



2\n+l 



2 (1 + X 2 ) 



2\n+l 



1 - (I + n)x 2 ' 



(28) 



where x = Rds/Rc- Considering R c of the order the 
curvature scalar today we have \i ~ 1. 

In the region of high density (i? R c ), we have 



2(2tt) 2 



d A kJk 



L eff 



(24) 



As we previously said, a regularization that do not re- 
spect the symmetries of the problem is incorrect [l8| , and 
produce the terms detailed beforehand. Hence a dimen- 
sional regularization of the energy density gives in the 
MS scheme 



4' c 



lim 

d-s-4 2(27r) d - 



jd-l 



kjk 



l eff 



"°eff 

64tt 2 



In 



( m eff \ 



(25) 



where fiQ is introduced to clean up the units. 



We shall hereafter would specialize to f(R) gravity. 
We should emphasis that the scalaron potential in gen- 
eral is a complicated one and certainly does not belong 
to the class of renormalizable theory. However, in the 
neighborhood of its minimum, the latter can be approxi- 
mated by a polynomial. Thus we can apply the quantum 



f(R) ~ R — plR c [l — (R/R c ) 



-2nl 



(29) 



It can easily be noticed from (|14[) that the minimum 
is R ~ p/Mpi as in General Relativity. Hence the gravita- 
tional sector is exactly equivalent to the standard frame 
work of General Relativity. The scalaron which settles 
at the minimum of the effective potential has small vari- 
ation around this point because of the space dependence 
of the density of matter (Ti"4")) . 

This translates to the chameleon field via its definition 
© and gives the minimum of the effective potential 



M, 



pl 



In /' 



f'(-?-) 



(30) 



Let us now consider the experimental bound that 
comes from the solar system tests of the equivalence prin- 
ciple (LLR). Using the thin-shell parameter [l4| for the 
Earth e t h we have 



eth 



6|/3|M pi $ e 



< 



x 10" 



(31) 
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where (</>oo,00) are respectively the minimum of the 
effective potential at infinity and inside the planet and 
$© the Newton potential for the Earth. 

Using the value $® ~ 7 x 1CP 10 , the previous bound 
translates into 4>oo/M p i < 10~ 15 , which after using 
Eq. <[30]) leads to 



fffSE.) _ i 



< 10 



-15 



(32) 



For the HSS model and with the density p c 

-3 



10 



-24 

g cm- and R c ~ H§, we have l0- 5 ( 2 ™ +1 ) < 10~ 15 tells 
us that n > 1 fill]. 

We will show that for this set of parameters the quan- 
tum stability is violated in the Starobinsky model as the 
mass of scalaron in the model grows fast with density and 
can easily cross the quantum bound. 



Quantum stability of f(R) gravity 

According to fl5j . the bound on m e // obtained using 
the 1-loop Coleman- Weinberg correction is given by 



m e ff(p) < 5.4 x 10" 



10 g cm 



-3 



1/3 



cV, (33) 



Also from Eq. ([15]) , we can express the scalar field mass 
in Starobinsky model as a function of the density p 



n+l 



m eff(p) 



v / 6pn{2n + 1) 



M 2 R C 



(34) 



where we assumed that the density is large enough 
compared to the cosmological density M 2 { R C ~ p c ~ 
10~ 29 g cm" 3 

From the previous discussion, we know that p ~ 1 and 
n > 1, which gives 



eff(p) 



3 x iQ-34 (£_\ eV . 



(35) 



At the cosmological density, p ~ p c ~ 10 _29 g/cm 3 , the 
quantum stability bound, 5 x 10 -13 eV, is larger than 
the scalar field mass m e ff = 3x 10~ 34 eV. However, the 
quantum bound oc p 1 ^ 3 while m e ff oc p n+1 (with n > 1). 
It is therefore clear that the m e f / will be excluded easily 
by this quantum stability bound at some high density. 
Indeed, the scalar mass m e ff is quickly excluded by the 
quantum stability bound at p rss i()- 87 ™/( 2 + 3 ") g/cm 3 . 

That means according to this bound, f(R) gravity in 
the point of view of chameleon theory is excluded in any 
typical dense medium, i.e., in the air (p air ~ 10 -3 g/cm 3 ). 



Extended Starobinsky model: Adding aR 2 term 

In high density regime the quantum curvature cor- 
rections to Einstein- Hilbert action become important. 
These corrections might provide a cut off to the scalaron 
mass. In what follows, we shall consider the extension of 
Starobinsky model by adding aR 2 term to its action, 



f(R) = R — pR c 



1-1 * 

Rc 



aR 2 , 



(36) 



This correction was briefly suggested in the original 
paper Q to avoid the problem of scalaron mass from be- 
coming too large and being inspired by the Starobinsky's 
original idea, it was introduced in [12j as a solution to 
the Frolov singularity problem(see also Refs.[l3| on the 
same theme) . As noticed above, the addition of this term 
does not change the position of the minimum of the field 
but it will provide a natural upper bound to the mass 
of the chameleon field. Thus the gravitational sector is 
unchanged compared to general Relativity, we have the 
same curvature scalar but the scalar sector is modified 
because of this additional term. The shape of the effec- 
tive potential is changed. 

In the regime of large densities, we have 



l eff 



6a (1 + 2ap/M 2 ) 



(37) 



We should emphasize that the scalaron is massless in 
the Einstein frame when the density diverge contrary to 
the Jordan frame where the effective mass goes to l/6a, 
this is because of the conformal factor and it will 
have no effect on the following discussion. In fact, we 
consider hereafter p ~ pi a b which gives in both frames 
ml s ~ l/6a. 

The classicality condition ([33| gives a lower bound a > 
6 x 10 3 cV~ 2 when according to the bound from Big 
Bang nucleosynthesis and CMB physics a <C 10 35 eV -2 
[l6| . But the tightest bound comes from the Eot-Wash 
experiments which implies that [2(| a < 4 x 10 4 cV~ 2 . 
We still have a range of viability of the model as soon 
as we add the quadratic curvature correction term in the 
action. 

Also it should be noticed that i? 2 -term is introduced 
here as a cure of late time cosmic dynamics. However, we 
know that R 2 can gives rise to inflation at early epochs. 
And if the model is used to describe an early acceler- 
ation phase we would have @, HJ a ~ 10~ 45 (a75o) 

cV~ 2 where N is the number of e-folds. This is certainly 
not be compatible with the quantum bound. We should, 
however, note that at high energies in the early universe, 
the quantum correction may be quite different than the 
one given by Coleman- Weinberg potential and the simple 
analysis presented here may not be valid in that regime. 
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II. CONCLUSION 

In this letter, we have investigated the issues of clas- 
sicality of scalarons in the Starobinsky of f(R) grav- 
ity model. The model is studied in the Einstein frame 
where we have General Relativity with a scalar field, 
the scalaron. At the densities studied, local analysis, 
we have shown that the curvature scalar is equivalent to 
the one in General Relativity. The scalar field appears 
as a Klein- Gordon massive field and can be quantized 
using the standard procedure. In this context, we have 
shown that the quantum bound on scalaron mass derived 
in [i~5| can be a tight constraint on f(R) dark energy 
models. Within the range of viability of the parameters 
of the model, that we derived, the quantum corrections 
are large for low densities. The scalaron masses increases 
very fast with medium density than the quantum bound 
on it. The mass of the scalaron is unbounded and can 
exceed the Planck mass at a reasonable value of the den- 
sity of the medium. Clearly, the model cannot be trusted 
in this case. We therefore need to introduce a cut off to 
the mass of scalaron as the expression of our ignorance 
and use the model below the cut off. 
In view of the aforesaid, we used the extended Starobin- 
sky model by adding quadratic curvature correction to 
the Lagrangian.We have shown that this term does not 
change the curvature scalar R ~ p/M^ but only effects 
the form of the potential V and therefore the mass of 
the scalaron. This term influences principally the scalar 
sector of the theory. It produces a natural bound on the 
mass of the scalaron that we constrained via the classi- 
cality condition and the Eot-Wash experiments. We have 
demonstrated that extended scenario is consistent with 
the quantum bound on the scalaron mass. 

It should also be noticed that in the regime of high den- 
sities or low densities but large scales, gravity is never 
weak and quantization of the scalaron might become 



complicated. As noticed by Starobinsky several years 
back, the quantum corrections in his model during in- 
flation are small[22j and hence there is no strong bound 
on the mass of scalaron from the classicality condition in 
this case. The analysis performed in this letter is done 
in a regime where the gravity is weak, density is low and 
scales are small. Hence the standard results of quantum 
field theory could be applied. Last but not least, we 
should clearly point out a subtle point of our analysis. 
To be fare, underlying argument regarding the quantum 
bound relies on the assumption of a semiclassical gravity 
in which chameleon field is quantum and gravitational 
field classical. But scalaron field has geometric origin 
which controls the curvature of space-time. To be sys- 
tematic, a full analysis should be performed. Thus the 
scalaron field theory, which is obviously not rcnormaliz- 
able, can be judiciously used below some ultraviolet cut 
off M uv . The quantum corrected effective Lagrangian 
contains a term proportional to rnj t fM^ v which might 
effect the analysis presented here[22|; in our opinion, the 
problem requires further investigation. 
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